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Do you like Pasta? Do you like Mathematics? 
Then you possibly will like 



Pasta by Design, George L. Legendre, Thames & Hudson, ISBN 9780500515808 
Deutsche Ausgabe: Pasta und Design, Springer Spektrum, ISBN 978-3-8274-2978-0 





Or do you prefer Gnocchi? 

Buon Appetito! 
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Liebe bUG-Mitglieder, 

In Osterreich haben die Sommerferien 
begonnen und da will ich gleich alien jenen, 
die schon oder demntichst Ferien haben, ei- 
nen schonen Sommer und erholsame Tage 
und Wochen wiinschen. 

Mit diesem bNL will ich vor allem „alte 
Schulden" einlosen. Wie so oft haben auch 
kurze Anfragen und Beitrtige ihre Eigendy- 
namik entwickelt und zu weiteren Untersu- 
chungen angespornt. Aus „ Building Tower s' 1 
habe ich die Aufforderung zu eigenen Ex- 
perimenten ebenso gerne aufgegriffen wie 
mich ..Fresnel Integrals ..." an friihere 
bNL-Beitrage zur Cornu Spirale und zu den 
„natiirlichen Gleichungen" erinnert haben. 
ba entwickelte sich gleich ein lebhafter 
Briefwechsel mit bavid Halprin ( Two good 
Turns), bavids interessante Laufbahn als 
Mathematiker verdient einen eigenen Bei- 
trag im nachsten bNL. 

Sie finden eine kurze Anmerkung von 
Erik van Lantschoot zu seinem Brusseler 
Tor in Dendermonde (bNL#89). Kiirzlich 
hat er eine weitere Ausarbeitung basierend 
auf seinen intensiven Nachforschungen in 
historischen Quellen geschickt. 

Christine Kova besuchte vor einigen 
Jahren einen DERIVE &. GeoGebra Work- 
shop “Fadengrafik". Spater schickte sie 
Bilder mit Mustern, die ihre Schulerinnen 
nicht nur am PC sondern auch mit Holz, Nti- 
geln und Faden erzeugt haben. Kurzlich er- 
hielt ich wieder Bilder, nun wagten sich die 
Schulerinnen bereits in den 3b-Raum. Ihre 
Skulpturen gegossen in Acrywiirfel werden 
sogar in einem Wiener Museum ausgestellt. 
Es ist schon, dass ein Mathematik- 
Workshop so nachhaltige Wirkung zeigen 
kann. 

Herbert Nieder verdanke ich den Hin- 
weis auf das „Pasta-Buch“. Obwohl ich kein 
besonderer Nudelfan bin, habe ich doch 
Appetit bekommen, die Agnolotti und Gnoc- 
chi zu modellieren. Im Internet konnen Sie 
noch einige Beispiele finden. 


bear bUG Members, 

Summer Holidays have started in Aus- 
tria. I'd like to wish all of you already being 
in holidays or waiting for them a wonderful 
summer and relaxing days and weeks. 

With this bNL I will settle “old debts". 
How so often short requests and mails de- 
velop their own dynamics and inspire for 
further investigations. I took seriously Ro- 
land's invitation in "Building Towers" to un- 
dertake own experiments and I was re- 
minded on earlier bNL-contributions on the 
Cornu Spiral and on "Natural Equations" 
reading “ Fresnel Integrals ..." . An animated 
communication with bavid Halprin followed 
( Two good Turns), bavid's rich career as 
mathematician deserves an own contribu- 
tion in the next bNL. 

You can find a short note from Erik van 
Lantschoot concerning his Brussels Gate in 
Dendermonde (bNL#89). Just recently he 
sent another paper based on his intense in- 
vestigations of historical resources. 

Christine Kova attended a bERIVE & 
GeoGebra workshop “Stitching Patterns" 
several years ago. Later she sent pictures 
with patterns produced by her students not 
only by means of mathematics but with 
wood, nails and threads. Three weeks ago I 
received some new pictures. Her students 
transferred the graphics into 3b-space. 
Their sculptures were cast in acryl cubes. 
Some of their works of (mathematical) art 
are exhibited in a Viennese museum. It is 
great that a mathematics-workshop can 
achieve such a sustained effect. 

We owe the note on the “Pasta-book" to 
Herbert Nieder from Hamburg. Although I 
am not a “noodles fan" at all, I got some ap- 
petite modelling Agnolotti and Gnocchi. You 
can find some more examples in the web. 

Viele GriiGe, Best regards 



Download all DTVL-Derive- and TI-f!les from 

http : / / www . aus tromath . at/ dug/ 
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The DERIVE-NEWSLETTER is the Bulle- 
tin of the DERIVE & CAS- 7/ User Group. 
It is published at least four times a year 
with a content of 40 pages minimum. The 
goals of the DNL are to enable the ex- 
change of experiences made with DERIVE, 
TI- CAS and other CAS as well to create a 
group to discuss the possibilities of new 
methodical and didactical manners in 
teaching mathematics. 


Editor: Mag. Josef Bohm 
D'Lust 1, A-3042 Wiirmla 
Austria 

Phone: ++43-(0)660 3136365 

e-mail: nojo.boehm@pgv.at 


I II 

Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. It 
must be said, though, that non-English ar- 
ticles will be warmly welcomed nonethe- 
less. Your contributions will be edited but 
not assessed. By submitting articles the au- 
thor gives his consent for reprinting it in 
the DNL. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE & CAS- 77 Newsletter will be. 


Next issue: September 20 1 3 


Preview: Contributions waiting to be published 

Some simulations of Random Experiments, J. Bohm, AUT, Lorenz Kopp, GER 
Wonderful World of Pedal Curves, J. Bohm, AUT 
Tools for 3D-Problems, P. Luke-Rosendahl, GER 
Hill-Encription, J. Bohm, AUT 

Simulating a Graphing Calculator in DERIVE, J. Bohm, AUT 

Do you know this? Cabri & CAS on PC and Handheld, W. Wegscheider, AUT 

An Interesting Problem with a Triangle, Steiner Point, P. Luke-Rosendahl, GER 

Graphics World, Currency Change, P. Charland, CAN 

Cubics, Quartics - Interesting features, T. Koller & J. Bohm, AUT 

Logos of Companies as an Inspiration for Math Teaching 

Exciting Surfaces in the FAZ / Pierre Charland 's Graphics Gallery 

BooleanPlots.mth, P. Schofield, UK 

Old traditional examples for a CAS - what's new? J. Bohm, AUT 

Truth Tables on the TI, M. R. Phillips, USA 

Where oh Where is It? (GPS with CAS), C. & P. Leinbach, USA 

Embroidery Patterns, H. Ludwig, GER 

Mandelbrot and Newton with DERIVE, Roman Hasek, CZK 

Tutorials for the NSpireCAS, G. Herweyers, BEL 

Some Projects with Students, R. Schroder, GER 

Dirac Algebra, Clifford Algebra, D. R. Lunsford, USA 

Treating Differential Equations (M. Beaudin, G. Piccard, Ch. Trottier), CAN 

A New Approach to Taylor Series, D. Oertel, GER 

Cesar Multiplication, G. Schodl, AUT 

Henon & Co; Find your very own Strange Attractor, J. Bohm, AUT 
Rational Hooks, J. Lechner, AUT 

Simulation of Dynamic Systems with various Tools, J. Bohm, AUT 
Belidor, E. v. Lantschoot 

and others 


Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wiirmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag. Josef Bohm 
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The third case x/O <y < Ox is still remaining. In this case we reduce [x,y] in both piles in such a way 
that the difference x - y equals the difference of a pair from v(m). So we have to substitute for 
else ??? the construct: 


IF 


^ x 

< y < cp ■ x a (v(m)) 

L <P 


n , 1 


- (v(m)) 

n , 2 


x - y, I" (v(m)) , (v(m)) 1, 0 

L n , 1 n , 2 J J 


The complete function reads now: 


#8: win(x, y, m) := VECTOR 


IF 


y > cp-x a (v(m)) 


= x, r x , (v(m)) 1, 

,1 L n,2j 


IF 


y < 


(v(m)) = y, rCv(m)) , yl, 

n , 2 L n ,1 J 


IF 


— < y < cp-x a (v(m)) - (v(m)) = x - y, 

V. <p n,l n,2 


rCvOtO) , Cv(nO) 

1. oil 

\ I 

L n.l 

n,2j )). 

J J 


#9: win (13 , 5, 10) 

#10: [0, 0, 0, 0, 0, 0, 0, [3, 5], 0, 0, 0, 0, 0, 0, 0, 0] 

We remove the zeros and make m dependent on x andy. Function wi n_move(x , y) is now ready: 
win_move(x, y) := (SELECT (w ± 0, w, win(x, y, MAX(x, y)))) 


#11 

#12 

#13 


win_move(13, 5) = [3, 5] 
win_move(30, 20) = 


[ 26 , 16 ] 


#13 needs 5.3 seconds calculation time. 


If our oppponent leaves piles of heights x and y we run wi n_move(x,y) in order to learn how to 
reduce the number of counters to gain the victory. In case of (13,5) we have to remove 10 counters 
from the higher pile and in case of (30,20) our strategy must be removing four counters from each 
pile. 

Two special cases are missing: in case x =y it is not necessary to find a DERIVE - answer because it is 
obvious that one will win with the next move. Case wi n_move(x , y) = [x , y] is only possible if the 
opponent starts the game and possibly knows the strategy. If so then we will loose - we can try our 
luck and make any move. 

#14: win_move(16, 26) = [16, 26] 

Agreeing upon y > x function wi n (x , y , m) (expression #8) can be simplified. 


Comment of the editor: I remember DNL#45 from 2002 when Richard Schorn wrote an article 
about Wythoff’s Nim. He closed his ACDC-contribution as follows: 
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The generation of the two sequences relies on some of the mentioned properties: 
Start with 1 as the x-value of the first safe pair. Add this to its position number to 
obtain 2 as the y-value. The x-value of the next pair is the smallest positive inte- 
ger not previously used, it is 3. The y-value is 5, the sum of 3 and the position 
number. The next x-value is 4, the smallest integer not yet used. 

My program for DERIVE is based on this algorithm in Martin Gardner’s book 
Penrose Tiles to Trapdoor Ciphers. (Richard Schorn) 

w(z, W := [[1, 1, 2]], x := 3, i := 2, M 
Prog 
Loop 

If i - 1 = z 
RETURN W 

W := APPEND(W, [[i, x, i + x]]) 

#1: M:=Mu{x,i + x} 

-i := -i + 1 
Loop 

If MEMBER?(x , H) 
x := x + 1 
exit 


See the TI-Nspire-version from 2013 which can 
of sets nor the useful MEMBER?-function: 




Now let’s proceed with Roland Schroder’s paper! 
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The author of this paper has - based on a small DERIVE program - discovered two more games with 
towers of counters which both are leading to the sequence <FLOOR(On) ne N>. Both games are played 
by only one player who intends to discover mathematical patterns in his results. Both games are 
starting with a sequence of piles of counters according to the figure given below. 



The sequence of 66 Countertowers ” or “Cointowers”can be continued to the left as far as one likes. The 
positions of the towers are numbered from right to left. When the game starts each tower consists of as 
many counters as its number says. For the first game we can assume an infinite number of piles. 

Every move consists of taking the last pile (first from right) and distributing its counters piece for 
piece on the next piles. We can describe the constellation at the beginning of the game as: 

... 7, 6, 5, 4, 3, 2, 1. Then we have: 

... 7, 6, 5, 4, 3, 3 after the first move and 

... 7, 6, 6, 5, 4 after the second move, etc. 

The next moves are left for the reader. (He/she might use the DERIVE program which will be 
presented at the end of the article.) Inspecting the piles to be distributed in the next move and counting 
its counters we find the sequence <1, 3, 4, 6, ...>. We can conjecture that here again the sequence of 
the first components of Wythoff s game appears. 

The second game deals with a finite number of towers. As the moves follow the same rules as given in 
the first game there will come a moment when no more piles are at the left end of the row to deposit 
counters. In this case we proceed by adding one-counter piles in left direction. See an example starting 
with 4 towers: 

4 , 3 , 2,1 
4 , 3,3 
1 , 5,4 
1 , 1 , 2,6 
1 , 1 , 1 , 2 , 2,3 

1 . 1 , 2 , 3 , 3 

1 . 2 , 3 , 4 


Having reached the reverse order all piles are moving to the left without changing the number 
sequence. It needed 6 moves to reach <1, 2, 3, 4> from <4, 3, 2, 1>. One can try starting with 5, 6, 7, 
... towers and count how many moves are necessary to reach the reverse constellation of the starting 
towers. Doing this with paper and pencil (or with the DERIVE-Qxogmm) you will get the following 
table: 


Number of towers at the beginning 

i 

2 

3 

4 

5 

6 

Number of moves to reach the mirror position 

i 

3 

4 

6 

8 

9 


This sequence in the 2 nd row is again very suspicious to form the sequence of the x- values from 
Wythoff s game! 
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We will simulate the second game using DERIVE. This simulation will give the results for the infinite 
form, too. Let’s start with an initial list (vector) representing the towers with heights from 1 to n. 


# 1 : init(n) := REVERSE (VECTOR (k, k, n)) 

h (z) removes the last element of any sequence of towers (the utmost right tower) and i (z) transforms 
this pile of counters in a list of single counters. 

#2: h(z) := REVERSE (REST (REVERSE (z))) 

#3: -i(z) := VECTORCL, m, FIRST (REVERSE (z))) 

(We use the list-oriented commands like REVERSE, FIRST and REST because DERIVE is programmed 
in LISP which is based on working internally with lists.) 

We rename the longer list as y and the shorter one as x: 

#4: [a := DIMENSION (h(z)), b := DIMENSION ( i (z)) ] 

#5: [y := IF(a > b, h(z), i(z)), x := IF(a < b, h(z), i(z))] 


The order of towers h(z) and i (z) are compared with respect to the number of towers. The greater 
number is y and the smalles one is x. y must be divided in two lists: in f (y) which consists of the same 
number of elements as x and g(y) containing the remaining elements. 

f(y) := VECTORCy , n, |a - b| + 1, DIMENSION Cy)) 

#6: n 

gCy) := VECTORCy , n, 1, |a - b|) 

#7: n 

Now we can add f (y) and x and append g(y) in order to get the next sequence of towers. 

#8: new_towersCz) := APPENDCgCy). f(y) + x) 

For better understanding the procedure we will demonstrate the process starting with a sequence of 
six piles: 

z := [1, 1, 1, 2, 2, 8] 

hCz) = [1, 1, 1, 2, 2] 

-iCz) = [1, 1, 1, 1, 1, 1, 1, 1] 

[a, b] = [5, 8] 
y = [1, 1, 1, 1, 1, 1, 1, 1] 

x = [1, 1, 1, 2, 2] 

fCy) = [1, 1, 1, 1, 1] 

gCy) = [l, l, l] 

fCy) + x = [2, 2, 2, 3, 3] 

APPENDCgCy), fCy) + x) = [1, 1, 1, 2, 2, 2, 3, 3] 

All steps are collected in one iterative procedure, starting with i ni t(n) : 

#9: towers (n) := ITERATES (new_towe rs (z) , z, init(rO) 


My first game starts with 5 towers: 
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towers(S) = [[5, 4, 3, 2, 1], [5, 4, 3, 3], [6, 5, 4], [1, 1, 7, 6], [1, 1, 1, 2, 2, 8], [1, 1, 
1, 2, 2, 2, 3, 3], [1, 1, 1, 2, 3, 3, 4], [1, 1, 2, 3, 4, 4], [1, 2, 3, 4, 5], [1, 2, 3, 4, 
5]] 

You can find our z from above as 5 th element and its follower. It needs 8 moves to come from 
[5, 4, 3, 2,1] to [1,2,3,4,5]. 

Question #1: Can we confirm our conjecture from above about the number of elements until reaching 
the reverse constellation? DERIVE’’ s iteration stops when two iterations remain the same. 

DIM(towers(10)) = 18 

#18: moves(n) := APPEND([1], VECTOR(DIM(towers(k)) - 2, k, 2, n)) 

moves(30) = [1, 3, 4, 6, 8, 9, 11, 12, 14, 16, 17, 19, 21, 22, 24, 25, 27, 29, 30, 32, 
33, 35, 37, 38, 40, 42, 43, 45, 46, 48] 


r 


VECTOR 


k-(l + 75) ) 


FLOOR 


, k, 1, 30 


[1, 3, 4, 6, 8, 9, 11, 12, 14, 16, 17, 19, 21, 22, 24, 25, 27, 29, 30, 32, 33, 35, 37, 
38, 40, 42, 43, 45, 46, 48] 


Please compare the sequences! 

Question #2: How can we derive the results of the infinite game from these results? 

From the above towers (5) -result we can read off the last numbers of the intermediate steps: 
[1,3, 4, 6, 8], Then we would need a 6 th and an 7 th pile of counters to proceed. The sequence needed is 
the increasing partial sequence of the last elements of the towers(n)-results. A short program can 
help: 

gamel(n, g, 1 := []) := 

Prog 

g := VECTOR(FIRST(REVERSE(v)), v, towers(n)) 

Loop 

1 := APPENDC1, [FIRST(g)]) 

If FIRST(g) > gr2 
RETURN 1 
g := REST(g) 

gamel(20) = [1, 3, 4, 6, 8, 9, 11, 12, 14, 16, 17, 19, 21, 22, 24, 25, 27, 29, 30, 32] 


Again we have the first components of Wythoff s game. 


The interested reader might have found out by himself that there might “sleep” some other exciting 
discoverings. It could be worthwhile to change the start conditions. What we showed here are only 
results which affirm the conjectures, there are no proofs. The proofs together with other rules can be 
found in Roland Schroder, “Der Goldene Schnitt in Tiirmen aus Spiels teinen”, in Mathematische 
Semesterberichte, Band 60, Heft 1, Springer 2013. 

http://scienceindex.com/stories/2251230/Der_Goldene_Schnitt_in_Trmen_aus_Spielsteinen.html 
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Another Game - Another Pattern? 


Roland’s last paragraph made me considering about “other exciting discoveries" . What to 
discover? So I had the idea to start with towers consisting of odd numbers of counters: 
[..., 7, 5, 3, 1] and playing the game from above: 

init2(n) := REVERSE (VECTOR (2 ■ k - 1, k, n)) 

-irrit2CL0) = [19, 17, 15, 13, 11, 9, 7, 5, 3, 1] 

towers2(n) := ITERATES ( new_towe rs (z ) , z, init2(n)) 
towers2(5) 


[[9, 7, 5, 3, 1], [9, 7, 5, 4], [1, 10, 8, 6], [1, 1, 1, 2, 11, 9], [1, 1, 
1, 1, 2, 2, 2, 3, 12], [1, 1, 1, 1, 2, 2, 2, 2, 3, 3, 3, 4], [1, 1, 1, 

1, 2, 2, 2, 3, 4, 4, 4], [1, 1, 1, 1, 2, 2, 3, 4, 5, 5], [1, 1, 1, 1, 3, 


[1, 2, 3, 3, 5, 5, 6] 


, [2, 3, 4, 


3, 4, 5, 6], [1, 1, 2, 2, 4, 4, 5, 6], 

4, 6, 6], [1, 3, 4, 5, 5, 7], [1, 1, 2, 4, 5, 6, 6], [2, 2, 3, 5, 6, 7], 
[1, 1, 3, 3, 4, 6, 7], [1, 2, 2, 4, 4, S, 7], [1, 2, 3, 3, 5, 5, 6]] 


We cannot find the reverse order of the initial status. But we notice a cycle of length 7 
starting with towers 1,2,3, 4, 5, 5, 6. 

We can define a function to calculate the length of the period: 

per_length(n) := DIM(towers2(n)) - P0SITI0N(FIRST(REVERSE(towers2(n))) , 
towers2(n)) 


per_length(5) = 7 

How long are the periods for games from 1 to 30 starting piles with odd numbers of counters: 
VECT0R([k, per_1ength(k)] , k, 1, 30)’ 

145 sec 


. 2 

3 

4 

5 6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

3 

4 

6 

7 1 

10 

11 

13 

14 

16 

17 

18 

20 

21 

23 

24 

25 

27 

28 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 ' 









30 

31 

33 

34 

35 

37 

38 

40 

41 

42 . 










5 piles lead to a period of length 7. I cannot detect any generating rule or pattern, can you? 
Next question: How many moves are necessary until having finished the first cycle? 
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moves2(n) := APPEND([[1, 1]], VECTOR([k, DIM(towers2(k)) - 2], k, 2, n))' 



' 1 

2 

3 

4 5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

IB 

moves2(20) = 

. 1 

5 

8 

12 16 

22 

23 

27 

31 

34 

146 

42 

209 

147 

199 


16 

17 

18 

19 

20 ‘ 

348 

86 

468 

279 

521 . 


Another function lists the length of the pre-periods: 

pre_per2(n) := VECTOR([k, P0SITI0N(FIRST(REVERSE(towers2(k))) , towers2(k)) 
- 1], k, 1, n)' 



1 2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

,per2(20) = 

1 

o 

UJ 

5 

7 

10 

22 

14 

17 

19 

21 

131 

26 

192 

128 


15 

16 

17 

18 

19 

20 ' 

179 

326 

63 

444 

253 

494 . 


I will check for the game starting with 6 piles form 1 , 3 11: 

towers2(6) = [[11, 9, 7, 5, 3, 1], [11, 9, 7, 5, 4], [12, 10, 8, 6], [1, 


1, 

1, 

13, 11, 9] , 

, [1, 

1, 1, 1, 2 

, 2, 

2, 14, 12], 

[1, 

1, 1, 

1, 2, 2, 

2, 

2, 

3, 

3, 3, 

15], 

[1, 1 

1, 1, 2, 

2, 2 

, 2, 3, 

3, 

3, 3 

’ > 4, 4 

, 4], [1, 

1, 

1, 

1, 

2, 2, 

2, 2, 

, 3, 3 

, 4, 4, 5, 

5], 

[1, 1. 1 

1 

, 2, 

2, 2, 

2, 4, 4, 

5, 

5, 

6], 

[1, 

1, 1, 

1, 2, 

2, 3, 3, 

5, 5, 

6, 6], 

[1, 

1, 

1, 1, 

2, 3, 4, 

4, 6 

6, 

7], 

[1, 

1, 1, 

2, 3, 

4, 5, 5, 

7, 7] 

, [1, 1, 

2, 

3, 

4, 5, 

6, 6, 8], 

[2, 

2, 

3, 

4, 5, 

6, 7 , 

, 7], 

[3, 3, 4, 

5, 6, 

7, 8], 

[1, 

1, 

4, 4, 

5, 6, 7, 

8], 

[1, 

2, 

2, 5 

, 5, ( 

5, 7, 

00 
i i 

i — i 

H 

l\J 

3, 3 

, 6, 6, 

7, 

8], 

[1, 2, 

3, 4, 4, 

7, 

7, 

8], 

[1, 

2, 3, 

4, 5, 

i — i 
00 

00 

i_n 

[1, 

2, 3, 4, 

5, 

6, 

6, 9], 

[1, 1, 2 

, 3, 

4, 

5, 

6, 7, 

7], [ 

:i, 2, 

3, 4, 5, 

6, 7, 

Efl, [1, 

2, 

3, 

4, 5, 

l 1 

1 1 

00 

<x> 



This is a strange result, indeed. The period has length 1 (compare with the list of the period 
lengths) and the pre-period consists of 22 elements. 


You are friendly invited to offer other ideas and games - and conjectures - and solutions - 
and proofs, Roland and Josef. 
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It was during Christmas holidays when my wife Noor and I made walk through 
winter wonder world. Suddenly the mobile phone interrupted the silence. At first I 
didn’t know who was calling. It was a colleague with whom I collaborated many 
years ago developing a curriculum for a special type of vocational schools. He 
was teacher at a Secondary school for forestry in Styria. 

After exchanging some personal small talk about ourselves (family, retirement, 
health, ...) he tried to explain his problem. I said, “Franz, send an email explain- 
ing your request and add one or two examples. I believe to understand you, but 
written text is quite better”. Coming home from our walk I turned on the notebook 
and found his mail: 

Franz Jauk’s Statistics Problem 


Dear Josef, 

... I have a small problem - in earlier times we produced our own programs ... 

The problem deals with finding a correlation coefficient. Given are a few pairs of values and their fre- 
quencies, e.g. (5,6) 25 times, (4,3) 17 times, etc. I don’t want to enter 25 times (5,6) and then 17 times 
(4,3), etc. I would like to enter 5,6,25 followed by 4,3,17, etc. That is it. 

Slope and intercept of the regression line are not so important. It would be great if I could represent 
the regression line but this is not a „must“ for me. 

Doing this in Excel I have to enter all pairs into the table - and this is boring. At the other hand the 
Excel 3D-representation looks fine: the value pairs in the xy-plane and the frequencies in z-direction. 




Please excuse that I am bothering you in your well-deserved retirement. 
Best regards 


Franz 
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Dear Franz, 

... Now to your problem: 

Nice to do it with other technologies (btw, it is not too boring working with Excel by copying 
down the respective cells). What you are expecting can easily be performed with Tl-Nspire. It 
can be done without any additional programming. 

I am quite sure that you can do it with the older Tl-handhelds (TI-92, V 200 and TI-84), too. 
But at the moment - I am on holidays - I don’t have these tools at my disposal. 

The screen shot below shows the result using one of your data sets. The regression line 
could easily be plotted together with the scatter diagram. What cannot be done is a 3D- 
representation. 



■ 

d 

id 


■ 

E 

F 

■ 

H 

♦ 




=TwoVar(a[],b[],e[]): CopyV 


=LinRegMx(a[],b[],e[] 


1 

1 

1 


20 

Title 

Statistiken mit zwei Varia... 

Title 

Lineare Regression ... 


2 

0 

0 

70 

X 

0*25 

RegEqn... 

m*x+b 


3 

0 

1 


5 

lx 

25. 

m 

0.733333 


4 

1 

0 


5 

lx 2 

25. 

b 

0*066667 


5 




SX := Sn-... 

0.435194 

r 2 

0.537778 


6 





ax := an... 

0.433013 

r 

0.733333 


7 





n 

100. 

Resid 

{0.2,-0.0666666666... 


8 




y 

0.25 




9 




sy 

25. 




10 




ly 2 

25. 




11 




sy := sn-... 

0.435194 




12 




ay := an... 

0.433013 




13 




Ixy 

20. 




14 




r 

0.733333 





r is the correlation coefficient and r 2 the coefficient of determination . 

You can also work with GeoGebra. But here you have to enter the data pairs (copying down 
the cells like with Excel). I did not find a way to work with triples (pairs together with their fre- 
quencies). 


franzj.ggb 
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The data set consists of 25 times (0,0), 20 times (1,1), 15 times (0,1) and 5 times (1,0). List 
Listel was created by a mouse click into the two spreadsheet columns. The correlation coef- 
ficient a = 0.57907 was calculated by CorrelationCoefficient(Listel) and the regression line 
b = -13x+24y=8 by FitLine(Listel) (in German: Korrelationskoeffizient(Listel) and Trend- 
linie(Listel)) 


It is not possible to have a 3D-representation of the 4 data piles. 


Let us turn to DERIVE. Two short programs produce the statistics and a 3D-plot as re- 
quested: 


reg_korr(md, mdext, d_, i, j, 11, 12, rl, r2) := 

Prog 

d_ := DIM(md) 
mdext := [] 
i := 1 
Loop 

If i > d_ exit 

j := 1 

Loop 

If j > mdj/i|3 exit 

mdext := APPENDfmdext, [[mdj,iil f mdj,ii2]]) 

j :+ 1 

i : + 1 

11 := FIT([x, a-x + b] , mdext) 

12 := FIT([x, a-x + b], S WAP_E L E ME NTS (mdext 1 , 1, 2)’) 

[rl := 0(11, x), r2 := 3(12, x)] 

[ "Regl: ", 11; ir Reg2: ", 12; r, r A 2: ", rl-r2; V: ", ^(rl-r2)] 

plot3d(md, mdext, pi, i, j) := 

Prog 

mdext := [] 
i := 1 
Loop 

If i > 4 exit 
pl := [] 
j := 1 
Loop 

If j > mdii|3 exit 

pl :=APPEND(pl, [[[md^iil, mdj,ii2,, j]]]) 
j :+ 1 

mdext := APPEND(mdext , [pl]) 
i :+ 1 
mdext 


See how they are working: 

(1 ) I enter the data as a matrix (named mm): 

1 1 55 

0 0 25 

0 1 15 

1 0 5 


(2) reg_korr(matri x) returns all interesting statistics: 


Regl: 

0. 5416666666- x + 0.375 

Reg2 : 

0. 619047619- x + 0.1666666666 

r A 2: 

0.3353174603 

r : 

0.5790660241 


You are invited to compare with the GeoGebra results from above. 
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(3) Apply plot3d(matrix) and receive a nice 3D-representation of the pillars. 

plot3d(mm) 



Back home some days later I could accomplish my proposals working with my good old Voy- 
age 200: 


The Stats/List Editor is the appropriate tool on the Voyage 200: 


r FIT Y F£t Y F]:t Y FHt ' 

|Tools|Plots|List|Calc 

F ST Y FbT y F?tY -i 

Distr|Tests|lnts| 1 

listl 

list2 

1 ist3 

1 ist4 

lists 

1 ist6 

1 

1 

55 




1 

G 

5 




G 

1 

15 




G 

G 






list3E4]=25 


HlilH FiliD AUTD PE 5/ ? 


FIt Y FEt T F5t 
T o ols|Flot s |L i st 


4: QuadReg... 
5: CubicReg... 
&: QuartReg... 
7: LnReg... 

S: ExpReg... 

0: PowerReg... 
fl: Logist83... 
0: Logistic... 
CIS in Reg 


f F Et Y FbT ¥ Y i 

□ istr|Tests|lnts| | 



lis 

1 

1 

G 

G 


i - ID i -+ - II n - 4- ir- - 1 - in ■; - + it-- 4- - i t .-4- - i 

t LinRtttax+M... V 

X List: 


1 t 

Y List: 

[lists 


Store RegEqn to: 

none* 


Freq: 

1 ist3 


Category List: 



Include Categories: 

O 


^CEnter=0K } £ESC=CFlNCEL^ 


t6 


list3E4]=25 

TVPE « [ENTER]=DK AND [ESC]=CflNCEL 


r FIT Y FEt Y FjT Y FHt Y F ST Y FbT Y Y 'i 

IToo 1 s Ip 1 ots II i st bale !□ l str It ests 1 1 nts 1 

list] 

f LinRt3(ox+t»... ^ 

i st6 

1 

1 

G 

G 

y=ax+b 

a =. 541666666667 

b =. 375GGQGQGQGQ 

r* =. 33531746Q317 

r =. 579Q66Q24144 

, CEnter=0K 5 ; 






list3[4]=25 

MAIN RAP AUTO PE 3/ ? 


Dear Franz, now it is your choice which tool to use. It was very interesting for me to compare 
the several possibilities treating your question. Now I have a question: Why do you need this 
tool in order to find the statistics for these special four data pairs? 


Josef 
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c) they are very susceptible for rounding errors when calculating greater values. So it is neces- 
sary evaluating them in exact mode only which causes extra long calculation times very soon. 

What to do? Although the tools given in the utility files are not very useful in my opinion, we can use 
DERIVE to calculate the integrals requested using another DERIVE feature. Luckily the Fresnel inte- 
grals can also be calculated applying the complex error function ERF(z). ERF can be calculated in two 
ways: as power series and as continued fraction. The first one converges very fast for small values and 
the second one converges fast for large values. And ... it is implemented as a standard function in 
DERIVE. This makes possible calculating Fresnel integrals and clothoids in real time with DERIVE 
without any extra programming efforts. 

I am sending a short file containing functions for Fresnel integrals and for clothoids. Maybe that one 
could produce an updated Fresnellntegrals .mth file. 

With best regards 
Piotr Trebisz 


Comment of the editor: I asked Piotr for possible applications of the Fresnel integrals (except 
the clothoid) and he answered: 

Unfortunately I cannot offer additional examples at the moment. It was the treatment of the clothoid 
which led me to these integrals. I know that they are used in physics for calculation of wave diffrac- 
tion, they are used in quantum mechanics and I found them in Richard Feyman’s path integrals. And 
as you wrote in your mail clothoids are used in road construction for forming smooth connections be- 
tween segments of different curvatures - and in particle accelerators. 

The Fresnel functions are the integrals of cos(a • t 2 ) and sin(a -t 2 ). a is an arbitrary real number. 


Let L the arc length of the clothoid and A its parameter (positive real number). Its curvature is propor- 
tional to Its asymptotic points are 


A-^ A-^fn 


and 


A-Jn A-^fn 


You can plot the clothoid expressions #5 or #8 of the following DERIVE file after introducing a slider 
bar for a in real time. 


With my best regards 
Piotr Trebisz 


What follows is Piotr’s DERIVE file. As I was very inspired by his mail and mth-file I did some 
research in my books and in the web. The list of my findings is given in the references at the 
end of this constribution, Josef 
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#1: [CaseMode := Sensitive, InputMode := Word] 

# 2 : 

#3: [a :€ Real (0, «>) , s :e Real] 

#4: 


#5: 


#6: 

#7: 

#8: 

#9: 

# 10 : 


*JjT • < 


KL0TH0IDE(a, s) := 

Jn-di C 

(1 + O-ERF 


(1 + l)-ERF 


s.(l - O 


+ (1 - O-ERF 


2- a 


s-(l + O 


2-a 


'N'N 


s-(l + O 


2-a 


+ (1 - O-ERF 


r s-(l - O 


Y\" 


2-a 


KL0TH0IDE_K0MPLEX(a , s) := KL0TH0IDE(a, s) • [1, -l] 
KLOTHOIDE_KOMPLEX(a , s) := 


r Jn-dL 

Jn- L-a N 

i 

-ERF 

' S.(l - i 

) N 

< 2 

i 

2 J 

< 2-a 

y 


[RE(KLOTHOIDE_KOMPLEX(a , s)) f IM(KLOTHOIDE_KOMPLEX(a , s))] 
lim [ RE (KLOTHOI DE_KOMPLEX(a , s)), IM (KLOTHOI DE_KOMPLEX(a , s))] = 


lim [RE (KLOTHOI DE_KOMPLEX(a , s)), IM (KLOTHOI DE_KOMPLEX(a , s))] 
s-»— ® 


*Jtt ■ a jJtt ■ 


2 2 _ 

■ a jJti ■ a 




d 

#12: — KLOTHOI DE (a, s) 

ds 



^ 2 ^ 


^ 2 ^ 


COS 

s 

, SIN 

s 


2 

2 





^ 2 ■ a j 


^ 2 ■ a j 



See now Piotr’s improved Fresnel integral functions together with their application defining 
the spiral in another way: 
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#14: 

#15: 

#16: 

#17: 

#18: 

#19: 


FRESNEL_C(z) := 


FRESNEI S(z) := 


(1 + O- ERF 


■Jtt-z-O- - i.) 


(1 + (.)■ ERF 


^7T-Z-(1 + i.) 


FRESNEL (z) := FRESNEI ,C(z) + i..FRESNEL_SCz) 


FRESNELCz) := 


■ ERF 


■Jtt-z-Q. - i.) 


a- [RE (FRESNEL (z)) , IM (FRESNEL fz)) ] 


lim a-[RECFRESNELCz)), IM (FRESNEL (z) ) ] 

Z-MO 




+ (1 - i_) ■ ERF 


v 

f 


+ (1 - i_) ■ ERF 


v 


J7T - Z ■ (1 + i.) 
2 

Jtt-z-CI - O 
2 


Jj 

jj 


a 


2 



Original DERIVE Utility: 


#24: FRESNEI SIN(2) = 



f 2 ^ 

SIN 

7T« t_ 

^ 2 ^ 



dt_ 


0 


#25: FRESNEI SIN(2) = 0.34341567S2 

0.14 sec calculation time 
Piotr's Fresnel SIN Integral: 

ERFO-CL - i.)) + ERF(>.(1 + i.)) L-CERFO-a + i.)) - ERFCjir-Cl - i.))) 
#26: FRESNEL_S(2) = + 


#27: FRESNEL_S(2) = 0.3434156783 


0.03 sec calculation time 
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It is interesting comparing the velocity plotting the spiral when using the DERIVE Utility func- 
tions - FR(z) - and Piotr’s functions - FRESNEL(z) with -2n <z< 2n. 

#28,: FR(z) := FRESNEL_C0S(z) + i ■ FRESNEL_SIN(z) 


#29: FR(z) 



r 2 1 



r 2 ^ 

COS 

7T- t_ 

■ -J 
+ 

1 

+- 1 

Tj 

SIN 

7T« t_ 




. 2 j 



. 2 j 


dt_ 


0 0 
#30: [RE(FR(z)), IM(FR(z))] 

#31: [RECFRESNELCz)), IM (FRESNEL (z) ) ] 

In Eric Weisstein’s Encyclopedia I found 
the plot of the slope of the Cornu Spiral: 


#32: 


FRESNEI S(z) 

FRESNEL_C(z) 


One of my favourite books dealing with 
differential geometry is Alfred Gray’s 
Modem Differential Geometry of Curves 
and Surfaces. 

Gray presents a generalization of this 
spiral with exchanged x- and /-coord- 
inates. cloth(l,a) is the clothoid 

S 


#33: clothfn, a, s) := a- 


#34: lim clothfl, a) 




r n + 1 1 



r n + 1 ^ 

SIN 

t 

+- 1 

Tj 

COS 

t 

^ r\ + 1 j 

^ n + 1 ^ 





dt 


Jtt ■ a Jtt ■ 


#33: clothfl, a) 


#36: cloth(2, a) 
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There is one chapter in Gray’s book treating Plotting plane curves with given curvature. Gray 
uses MATHEMATICA and offers a special tool to find the natural or intrinsic equation 1 * 1 of 
curves solving a system of DEs. 


I cannot transfer this tool directly to DERIVE but I can reproduce a nice Clothoid Flower. 


#37: cloth2(n, a, s) := a- 


s 



s 



m 

f n + 1 1 


’ 

r n + 1 1 


COS 

t 

dt, 

SIN 

t 

dt 




^ n + 1 ^ 



^ n + 1 j 


0 



0 





r 

' cos(cO 

SIN(cO ' 

2-7T 7 T ^ 

#3S: VECTOR 

cloth2(l, a) ■ 



i GE| 0, , 



. - SIN(cO 

cos(cO . 

3 3 J 



And here is another one (which is not shown in Gray’s book): 


#39: VECTOR 

r 

cloth2(2, 1.5)« 

' cos (a) 

SI NO) ' 

7T ^ 

, ffl, 0, 2 ■ 7T , 



. - SI IMCa) 

COSO) . 

4 J 



w Find more about the Natural or Intrinsic Equation in DNLs#29 and #30 (David Halprin’s 
articles). See also the rich list of references at the end of this contribution, Josef 
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I was so much inspired by some of the references - especially by [4] and [11] that I tried to 
find an easy to use DERIVE implementation for plotting curves with given functions of their 
curvature dependent on their curve length. Josef 
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Natural (Intrinsic) Equations of Curves 

Josef Bohm, Wurmla, Austria 


Hans Walser presents in his paper [11] a Turtle-Graph method using MuPad to find curves with given 
initial point, given start direction and a function k(s) describing the change of direction of the tangent 
dependent on the arc length s. k(s) is the curvature, its reciprocal value is the radius of the osculation 
circle. 

Walser starts with an introductory example - plotting a pentagon: 


poly_turt(n, 1, a, start, xo, yo, dir, xn, yn, pts, i) 
Prog 

xo := startj.1 
yo := startj.2 
dir := startj.3 
pts := [[xo, yo]] 

-i := 1 
Loop 

#1: If i > ri 

RETURN pts 

[xri := xo + l-COS(dir), yn := yo + l-SIN(dir)] 
pts := APPEND(pts, [[xn, yn]]) 
dir := dir + a 
xo := xn 
yo := yn 

i :+ 1 


# 2 : 


po]y_turt 


5, 2, 


2 «rr 

, [ 0 , 0 , 0 ] 

5 


We perform n moves (turning n times to the left by angle a and proceeding the distance / in the now 
new given direction, start gives the x- and y-coordinates of the starting point and as third component 
the initial direction. Repeating turning and proceeding 5 times results in a pentagon 



We start at any position and move to position s where we have to change our direction according a 
function of s. Performing a numerical approximation we can do this only executing very small steps 
c/s. It turns out to be something like a numerical approximation of a differential equation: 
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nat_equ(K, 1, ds, start, n, xo, yo, dir, xn, yn, pts, 
Prog 

xo := start^l 
yo := startj,2 
dir := startj,3 
pts := [[xo, yo]] 
n := ]/ds 

i := 1 

Loop 

#3: If i > n 

RETURN pts 

xn := xo + ds-COSCdir) 

yn := yo + ds'SIN(dir) 

pts := APPENDCpts, [[xn, yn]]) 

dir := dir + ds-SUBST(K, s, i-ds) 

xo := xn 

yo := yn 

i :+ 1 


k(s) is the curvature function of the arc length, / is the total distance on the curve, ds is the step 
width and start are the initial values as explained above. 


Take the first example: What is the curve with 
curvature = 2 in all of its points? I am quite sure 
that you know this: 

#4: nat_equ(2 , 2, 0.001, [0, 0, 0]) 


#5: x + 


1 V 

2 j 


1 

4 


We start at (0,0) looking in x-direction and 
change the direction after each step in the same 
constant way, i.e. by the value 2. This gives a 
circle with radius 0.5. You see the result of 
nat_equ and superimposed the respective circle 
(in grey colour). My approximation fits pretty 
well! I am satisfied. 



We started with Piotr's Clothoid. The curvature is 
proportinal to the distance made: 

#6: nat_equ(s , 6, 0.01, [0, 0, 0]) 

One could add any parameter a and then plot the 
curve like 

nat_equ(0 . 5 ■ s , 6, 0.01, [0, 0, 0]) 



0.2 0.4 0.6 0.S 


1 


1.2 1.4 
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nat_equ() invites to reproduce some of Gray's examples and to experiment with own functions k(s) 
Our reward are very strange looking figures as can be seen on the next page. 


2 

nat_equ(s , 6, 0.01, [0, 0, 0]) 



Do you remember Piotr's explanation that the 
Clothoid is the equivalent to the logarithmic 
spiral? 


nat_equ(s + SIN(s), 10, 0.005, [0, 0, 0]) 
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2 

nat_equ(s-SIN(s) , 20, 0,01, [0, 0, 0]) 


2 

nat_equ(s-SIN(s) , -20, -0,01, [0, 0, 0]) 



And this is my final natural equation curve: 


nat_equ(2.s.COS(s).SIN(s), 20, 0,005, [1, 1, 0]) 



Many thanks to Piotr Trebisz for his mail about improving the Fresnel Integrals, Josef. 
Can we do this with Tl-Nspire, too? 
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Natural Equations of Curves & Tl-Nspire 

It is not too difficult to transfer the D£/?/\/£-treatment on the Tl-Nspire environment. What 
we cannot do is defining the Fresnel integrals supported by the error function because it is 
not implemented there. 


fresnel_sin\z}\= 

W 

z 

sin 

0 

2 

JT U 
2 . 

d u \fresnel_cos\zj : = 

■ / 

2 

COS 

0 

|f 2' 
ic • u 

\ 2 , 

Done 

du 

fresnel_sin{2) 






0.3434156784 

do th !/j : = \fres n els in[ t) jresnel_cos\t) ] 



Done 

dothi lOj 




[0.4681699786 0.4998986942] 


dothi 20) 




[0.4840845359 0.4999873350] 


doth[ 50 J 




[0.4936338026 0.4999991894] 


dothi 1 OO) 




[0.4968169005 0.4999998986. 


E 


7/99 


We cannot find the asymptotic center of the spiral by calculating the limit. But we see that 
the curve approaches (0.5, 0.5). Plotting the clothoid needs some time but it works - on the 
PC. I didn't try on the handheld. 



x2l ')=fresnel_sin(?) 
y 2 ( ) =fr e snel_c o s ( ?) 


0.1 


0.1 
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The program to plot the natural or intrinsic equation is quite identical to the DERIVE 
procedure. The lists of the coordinates are stored as lx and ly. I stored them under new 
names in order to produce a plot of three curves on same axes. This works in a very 
reasonable time. 


nat_equ[ 2, 2,0, 01,(0, 0,0 }) 

lx and ly are ready to be plotted! 

Done 

Ix_circ\=fx 

[o 0.01 0.019998 0.02999 0.039972 0.04994 0.0P 


A 


nat_equ 

De fine nat_e qu( kapp a l 1 , ds f sta rt ) = 

Prgm 

Loca 1 n } xoyO, di r, xnyn, i 

xo\ = {stc i r/[ 1 ] } yo : = { sta r/[ 2 ] } : di r : =sta r/[ 3 ] 

lx\={xo}\fy\={yo} 


I 


fy_circ\=fy 

[o 0. 0.0002 0.0006 0.0012 0.001999 0.002997 ► 
6 , 0 . 005 ,( 0 , 0 , 0 } ) 

lx and ly are ready to be plotted! 

Done 

Ix_doth\=lx 

[0 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.0> 


fy_cIoth\=fy 

[o 0. 1.25E-7 5.E-7 0.000001 0.000002 0.000004 ► 


i 1 

7 i n 

nat_equ\ — ^10,0.01/ 

0,0- 

\s 

2 \j 


n \= — 
ds 

For i,l,n 

xn\=xo+ds ■ cos [dir] 
yn\=yo+ds- sin(t/;r) 

Ix\ =augment(/x, { x/i } ) 

£y:=augment(/y, \yn } ) 
dir\=dir+ds- lim _ [kappa) 
s^i- ds u 
xo : =xn yo : =yn 
EndFor 

Disp "lx and ly are ready to be plotted! " 
EndPrgm 


lx and ly are ready to be plotted! 


Done 


Ix_ex :=!x: Iy_ex : = ly 

r n n n i n n i 40 n n 1 c 1 n 7 n n 1 c n n n n q q -■ n > v 

8/99 I 


1/13 
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A Student’s Problem with Homework Files 


... My name is Michael Klamecki. My teacher, Prof. Haller recommended to contact you because of 
two nasty DERIVE problems which I came across just recently. 

My first problem: 

I am able to open DERIVE on my notebook only as administrator. Opening it just the „normal way“ 
and trying to solve any equation it does not work at all: the window becomes blue and I am unable to 
perform any mouse click. 

My second problem: 

It happened when I loaded a previously saved DERIVE file. In a small pop-up window an error mes- 
sage appeared reporting that there is an error „parsing expression #7“. After confirming the Ok-button 
another error was reported and finally the file was opened. 


It would be great if you could give any advice. 

with best regards 
Michael 

DNL: 

Dear Michael, 

What concerns your “administrator-problem”, I asked Gunter Schodl, the DERIVE-expert. 
This is his answer: 


This is a problem of the administrator rights. I know this problem it appears sometimes when 
students install DERIVE on their laptops. My standard advice is: 

• Uninstall DERIVE completely (including all remains) 

• Install again (right click ... perform as administrator). 

You can also try to change the rights in the DERIVE folder for “All Users”. 

If this does not work: 

Look for another laptop where DERIVE works as you would like, remove DERIVE from your 
device and then copy it back (using an USB-stick) from the other laptop. This is how we 
could resolve this nasty problem. 

I don’t know the reason for this problem. Probably it is connected with the security settings of 
some computers. I didn’t investigate the reason, I have solved the problem in most cases. 

Hope this helps, many thanks to Gunter. 

I cannot give any advice for your second problem (the parsing error) without knowing the re- 
spective file. So please send the file. 

Regards 

Josef 
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Michael came back to me: 


... I uninstalled DERIVE manually and removed all respective files (except my homework files). Then 
I reinstalled DERIVE but the same problem occurred. Ok, I will work with DERIVE as administrator 
in the future. 

I attach two screenshots and the DERIVE file causing the “parsing error” in expressions #7 and #8. 
I noticed that the fractions in #8 are printed in an unusual way and the AND-operator is written as 
AND (A) . 



1) Der geradl "irrige Tell der Autobahn, in GeoGebra (PQ) 

#3: P := [0, 975] 

#4: Q := [1000, 0] 

#5 : g(x) := k(x) + d 

#6: [g(0) = 975, g(1000) = 0] 

#7: 50LVE([g(0)=975 , g(1000)=0] , [d , k] , Real) 

#3: [d=975 AND k=- 39/40] 

Best regards 
Michael 


DNL: 

Dear Michael, 

the error was easy to find - if one knows how to detect it. 

Look at expression #5: You use variable k for the slope in the linear function g(x). This 
should read g(x) := k x + d. But #5: g(x) := k(x) + d. Something must have happened. 

Copy the “wrong” expressions #7 and #8 into the Entry Line, then you will notice that they are 
under quotes, i.e. they are not mathematical expressions but have changed to strings. This is 
the reason that the fractions are looking strange. See my screenshot below: 

Why that? Go to expression #31 ! What have you done there? 


ITERATES Again! 
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^ X |"SOLVE([g(0)=975 , g(1000)=0] , [d, k] , Real) 1 ' 

#1: CaseMode := Sensitive 

#2: InputMode := Word 

1]) Der geradlinige Teil der Autobahn, in GeoGebra (PQ) 


#3: 

P := [0, 975] 


#4: 

Q := [1000, 0] 


#5: 

g(x) := k(x) + d 


#6: 

[g(0) = 975, g(lOOO) = 

0] 

#7: 

■SOLVE ([q(G)=975 , g(1000; 

>=0] 


2 

#31: k(x) := a>x + b-x + c 

#32: [k(402) = 626, k(571) = 695, k'(571) = h'(571)] 

#33: SOLVE ([k (402) = 626, k(571) = 695, k'(571) = h'(571)], [a, b, c], Real) 


You defined a quadratic function k(x). This is ok at this position and you will not face any 
problems saving the file. 

Reloading this file DERIVE finds k(x) among the stored functions. Performing then all com- 
mands from the very beginning DERIVE substitutes k x by k(x) which leads immediately to 
problems when solving the system in #7. 

I renamed k(x) as kl(x), stored the file and everything works. So take care in the future nam- 
ing variables and functions. I recommend using more-letters-function names, like curve(x), 
mo_way(x), kul(x), ku2(x), etc. 

I wish much success in mathematics and much fun with CAS. 

Best regards to you and to my colleague Wilhelm Haller, 

Josef 

ITERATES again! (from Michel Beaudin) 

Dear Josef, as you showed in the panel discussion in Estonia, the « iterates » function of Derive works 
well when exported to Nspire CAS (in fact, the first export was at the Bonn conference as you wrote 
in DNL). But it seems to be less general than the one Derive has. 

For example, in Derive, iterates(diff(f,x), f, x A 10, 3) returns [x A 10, 10x A 9, 90x A 8, 720x A 7]. In Nspire 
CAS, the SAME command gives [x A 10, 0, 0, 0]. It seems that Nspire CAS computes diff(f, x) BEFORE 
replacing f by x A 10. Do you have an explanation? 

The reason for my question is the following : I have defined many functions in Nspire CAS and, in 
some cases, I was unable to use the iterates function. So I have used recursion which becomes slow. 
With iterates, it would be so fast! But I am not a good programmer compared to you. So, maybe 
you can help me! 

Best regards and hope to see you and Noor soon, 

Michel 



P 32 


ITERATES Again! 


D-N-L#90 


Dear Michel, 

this seems to be really difficult. 

I tried several approaches. The problem is - as you mentioned - the immediate evaluation of 
the differential operator. When we call dif(f,x) Nspire evaluates to zero. 

We see the bold “d” which cannot be transferred into the program as unevaluated operator. 

I have a - very weak - solution writing the differential quotient into the program. 

Maybe that experts like Philippe Fortin will provide a more general solution. 


newton _its [f,v,vO,n):=iterates 


fw 


,v,vO,n 


Done 1 iterates 


0 . 2 ,^ 


{ 0.2, 5. 1,2. 74608, 1.73719,1. 44424, 1.4145^ 


ne\vton_its\u -2000. 


i,w,2.,2C|l[2l] 


4.573050519 


2000 ' 


, 0.2 


4.573050519 


© The Tribonacci Sequence iterative 


trib_seq («):=( itemtes f { w[ 2] , w[ 3], n[ 1 ]+m{ 2]+h{ 3] } , w,{ 1 , 1 , 2} ,n~ l.)) T [ 1 ] 


Done 


nib seq f 1 5j 




[l 1 

2 4 7 13 24 44 81 149 274 504 

927 

1705 3136] 


vO, 

/ 

n=0 

newton\f, v, vO, n) : = 

newton[f,v,vO,n-\)~ lim 

,n> 0 

v newton(f,v,vO,n- 1) 

— W 




dv j 



Done 


newton[x ^ -2000, v, 2 . , 2o) 


4.573050519 


© The Hennite polynomials 


henns(n):=iterates ({>*{2], 2- .v n{2]-2- (>v{ 3]-El ) - w{ 1 ] ,>v[ 3]-El } , w,{ 1,2- .v,o},n) 


Done 


hetms(4j 


4-.v 2 -2 
4-.V (2- A' 2 — 3) 
16 a 4 -48 v 2 +12 
x ■ I A 


4- x -2 
4-.v (2- jv 2 — 3) 

16-.v 4 -48 a 2 +12 8-.v l4-.v 4 -20-.v 2 +15) 4. 


Define iterates( vect, var, stall , numb ) = 

Func 

Local inc, in , iterlist, numvar, tvect 

If left(string(var), l)=" { " or leftl|string(s/ar/),l)=" { " : Goto vector 
iterlist = {start} 

For inc, 1, numb 

start: = \ect\ vat-start 

iterlist: =augment(iterlist, { start } ) 

EndFor 

Goto end 

Lbl vector 

iterlist:=start 

For inc, l, numb 

statt = vect \ vat -start 

iterlist: =augmenx(iterlist , start) 

EndFor 

list^matj iterlist, dim[startjj 
, Lbl end 


6/1 6j 


1/14 


iterates_f 


Define LibPub iterates_f(^?,«)= 
Func 

Local inc, iter list 
iterlist:={f0 } 

For inc,l,n 

fo-=j-[fo) 

dx 

iterlist: =augment(iter/B/, {/ft})| 

EndFor 

iterlist 

EndFunc 


My iterates_f does what you asked. 


iterates 

(-v 10 j 

1 

, 10 ,10- ,r 9 ,90- ,r 8 ,720- ,r 7 ,5040- v 6 J 

iiemies _/‘(sin(2\ 
■{ sin( 2- 

v),5) 

■ x),2- cos(2-.f),-4- sin(2- _v) , _ 8 - cos(2- v) ,16- sin(2-_r),32- cos(2-.v) } 


Regards as ever, Josef 


From Michel 

Thank you dear Josef. Your solution is not so weak and can be useful! 

This example shows (again) how Derive was special... And this kind of example (and others) 
could make a good presentation for the Derive/Nspire session in Malaga. Is there any chance you at- 
tend ACA 2013? If not, let us promise to have a good beer together in Austria in 2014! 

Michel 
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Christoph Kuderli's problem with Error Trapping (TI-NspireCAS) 


In a TI-NspireCAS program I tried to trap the error „Division by 0“ by the construct 
Try... Else... EndTry. This does not work because the Else-block will not be executed. 

I informed in the online manual and didn’t find any advice and I had no success with a Google re- 
search. Trapping other errors - e.g. dimension error - is possible applying Try ... Else ... EndTry. 

Can you please forward this deficiency to TI? I don’t find any possibility to do this by myself. 


probQ§ 


0 


w= 1 
s= 0 

17 

pmax= — 
23 


Done 


probO 

Define probo()= 

Prgm 

P := ® 

p2 := 0 

For i, 0,12 
For j, 0,12 

If (;# 0 or jto) and (i# 12 or j#12) Then 


13/18 


p2\-- 


i 12 -i 

1 

i+j 24 -i-j 


Endlf 

If p2>p Then 

p:=p2 

w :=i 

Endlf| 

EndFor 
EndFor 
Disp "w=" J w 
Disp "s =" ,s 
Disp "pmax^,p 
1/99 |EndFrgm 




See the program without (above) and with the - non working - error trapping routine (below). 



prob 


Define prob()= 

Prgm 

Local p,p2,ij,w,s 

p:=0 

For 7,0,12 
For 7,0,12 
Try 

7 12-7 


20 / 20 | 


p2:= 

Else 


i+j 24-i-j 
2 


p2:=- 

4 

ClrErr 
EndTry 
If p2>p Then 

p:=p2 

w:=7 

s:=j 
Endlf 
EndFor 
EndFor 
Disp "w=",w 
Disp "s=",s 
Disp M pmax=",^ 
EndPrgm 


s 
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I sent a mail to our TI-NspireCAS expert, Philippe Fortin: 

Dear Philippe, 

How are you? Here I am again with a Nspire-question (provided by a DUG-member): 
Is it possible to catch a “Division by Zero” error using the Try - EndTry procedure? 

So many errors can be handled by Try - EndTry. How to do with this one? 

I attach the respective program. 

It is no problem to circumvent the problem using an IF-contruction. 

Best regards 
Josef 


This was Philippe’s answer: 

I hope everything is fine for you ! 

In fact, computations as 1/0 or 0/0 doesn't cause an error... we do get a result, equal to 
undef. There must be some reasons for this choice (but I am not sure of which ones !). 

The question is then : is it possible to check that p2 contains the special value undef 

It is not possible to use 

if p2=undef then... 

If p2 is actually equal to undef, the result of this test will be undef = undef, which is not true 
or false 

It would be possible to use something like 

If when(p2 = undef, false, false, true) then... 

Another solution is to test the equivalent string value... 

If string(p2) = "undef" then... 

Last, since the expected result should be a number, it is possible here to test 
If getType(p2) = "EXPR" then... 

But this would not work if p2 could possibly contains values as sqrt(2) which also have the 
type EXPR... 

None of these solutions is perfect (all are using odd constructions... and the last one is not 
correct for all situations). It would be better to have a function directly testing for a undef 
value. 

Of course, a Try ... EndTry including the test statement if p2>p then... would also be a possi- 
ble approach... 

Best regards, 

Philippe 
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Dear Philippe, 

many thanks for your extended answer. 

Like you I tried many approaches ... but not the last one. It works. 

I will publish your answer + solution of the problem in the next DNL. 

It would be interesting to learn why “Division by 0” is a “special error” which needs an extra 
treatment. 

Best regards and thank you once more, 

Josef 


And here is Philippe’s final comment: 

Dear Joseph, 

Just a thought concerning the use of undef (which is not an Error). It could be related to 

limit computation... 

When entering 

lim(l/x,x,0) 

it is rather a good point to get undef instead of an error message... 
by the way, you probably noticed that 

1/0 returns undef 
1/0 A 2 returns infinity 

This is also consistent with limit computation... 

(this could be included in the answer you plan to publish !) 

Best regards, 

Philippe 


An additional note on Brussels Gate sent by Erik van Lantschoot: 


The Brussels Gate does not end fascinating me; at first I wanted to learn more about the Dutch 
engineers who built all three gates (Brussels, Mechelen and Gent - the last one has been destroyed) in 
1822. In the meanwhile I found out that the head of the engineers was a certain Cornelis Alewyn. 
Then I wanted to lear from where C. A. has got his experience. Who knows history of tec hn ics knows 
that one has to do some research in France because the succeeders of Vauban started applying 
algebraic methods for calculation and planning military buildings. And I was lucky! ! The inventor of 
the draw bridge with a counterweight on a curved track is Bernard Forest de Belidor (1697-1761) who 
gave the equation of the track as a formula. Applications can be found in the United States (two 
double line road bridges) and in four or five French towns and in Switzerland. I will prepare a 
respective paper investigatig the difference between this curve and the parabola which I calculated in 
Dendermonde. 


Best regards 
EvL 
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Two good Turns 

David Halprin, Australia 

This was written in 1992 and made use of Mathcad 2.0. forgive this transgression. 

There are many traps for the unwary when attempting to find a functional relationship between 
two variables so as to represent a plane curve in some coordinate system. 

Since the Cartesian framework is the usual default coordinate system let us examine a common 
pitfall. Many curves have a symmetry about one axis, (at least), and if one requires the simplest equa- 
tion to represent the curve, then perhaps one arc of the curve will be preferable, since the symmetrical 
arc will follow with ease. 

A trivial example is the circle. The simplest equation in Cartesian Coordinates is for a semi- 
circle, where we have y = f(x), whereas for the whole circle we have to combine it with y = -f(x) to 
have the expression [y - f(x)][y + f(x)] = 0 which gives us the usual quadratic form. One should note 
well that:- 

1) Such a quadratic form allows 2 values for y for any value of x as it does also allow 2 values for 

x for any value of y. 

2,2 2 

x = c . 

2) The circle is a composite curve, when expressed as y = f(x). If we had chosen a Polar Coordi- 
nate System then we would have a much simpler functional relationship between two variables 

and the circle is a simple curve of one arc, r = c 

One can manipulate symbols in such a way that two different curves, or two arches, or branches, 
of related curves, (reflected, rotated and/or translated), can be combined in one equation to make it 
look like one curve, when it may be a simpler exercise, and much simpler equation, to analyse it as a 
composite curve. For instance any curve which lies above the X axis may be conjoined with its reflec- 
tion beneath the X axis and the resultant equation may be encompassing the two parts as though they 
are two branches of the same curve, but the equation will be of a higher degree, and it is far easier to 
look at each branch separately. 

In the problem of the Railway Transition Curve, one can fall into the trap of assuming that the 
entire curve is best represented: - 

1) in Cartesian Coordinates and 

2) in one equation. 

However, staying with Cartesian Coordinates, and looking for only one of the symmetrical arcs, 
one would have obtained a far simpler solution, even more significant than in the case of the circle ex- 
ample. 

One necessary quality of a point of inflection is that it has no curvature at the point itself, but on 
either side of that point the curvature takes opposite directions, meaning that a concave arc changes to 
a convex arc. 

The force on the flange of the train's wheels has to be constant as it rounds the bend. The only 
way to state this with exactitude, is to talk about the relationship of the curvature to the arc length, 
namely: 
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A Railway Transition Curve, which is also used by road builders (road engineers), who want it 
to be such that the radius of curvature does not change discontinuously, since this would involve dis- 
continuities in the centrifugal (centripetal) forces between the rails and wheel flanges of the trains, or 
the tyres and road for cars and trucks. So what is required is that the radius of curvature is a continu- 
ous function of the arc length. 

The simplest and most natural functional relationship is that curvature (the reciprocal of the ra- 
dius of curvature) is proportional to the arc length. 

1 

K = — OC S. 

P 

This equation uses two intrinsic variables, which are not to an external framework, as are the ex- 
trinsic Cartesian (x,y) or the extrinsic Polar (r,0). There are many advantages to be gained from famili- 
arity with the intrinsic variables, of which there are three, from which we may choose any pair. 

Radius of Curvature, rho (p) 

Arc length, (s) 

Tangential Angle, phi (cp), where 

ds dx dy dx dy 

p = — , — = coscp, — = smcp, — = p • coscp, — = p • sincp. 

dy ds ds dy dy 

This curve has been studied for some time under various and has also proven to have been use- 
ful in optics. It is known as Klothoid, Clothoid, Cornu's Spiral, Euler's Spiral and Railway Transition 
Curve. 

In addition to the use of a pair of variables in a functional relationship to define a plane curve, 

such an equation can be split into two equations by the use of a parameter, which is often referred to as 

a dummy variable. When one eliminates the parameter from the pair of equations, one is back to the 
original equation in two variables. However, if one looks for a geometrical image for the parameter, it 
is often an interesting finding. Commonly, when transforming coordinates from one system to another, 
one is left with two parametric equations and the parameter is one of the variables associated with the 
coordinate system, from which one is transforming. 

e.g. Going from polar to Cartesian, the polar angle, 0 is commonly the parameter, and it is a very 
convenient one, especially when instructing a plotting program to draw the curve, since one tells it to 
range 0 from 0 to 2n and the curve is drawn with great facility. 

However, there are serious pitfalls to acceptance without inspection of any and every such pa- 
rameter, that may appear when transforming, as described, depending on the particular curve and the 
geometrical meaning of the variable, pertaining to that particular curve. 

e.g. When transforming from Intrinsic to Cartesian Coordinates, there are s, cp and p. Usually 
there is no problem with s and cp, since with most curves one can range their values fairly well, but the 
radius of curvature may not range, and hence one may obtain, at most, a small arc of the curve 
and, at worst, an impossible representation of the original curve. 

There are certain differential equations, involving elliptic integrals, called Fresnel Integrals, which 
have no closed form of solution, but which can be approximated with a strange looking series. These 
integrals can be used as parametric pair of equations to define the Clothoid. 
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In fact, since it is the clothoid and can be specified very in the two variables, s (arclength) and p 
(radius of curvature), it is only a matter of numerical approximation, which will find the point of in- 
flection. Further, the curve is compounded from two specially selected arcs of the clothoid, in the first 
and second quadrants respectively, and the omitted parts of the curve from the diagram are useless to 
railway engineers, since they are spirals about asymptotic points. The curve, as you will see below, is 
made up of two spirals, which are interconnected through a point of inflection, so the arc chosen is 
anywhere from the point of inflection up to as much of the spiral as it suits the engineers to use. 


Family of Railway Transition Curves in Intrinsic Coordinates 


1. 


1 ^9 

— = k = — 1 - = 2 Ks 

p ds 


Cesaro Intrinsic Equation 


2. cp = K • s 2 + L Whewell Intrinsic Equation 

(Let 5 = 0 at cp = 0 :.L- 0) 


9 


3. s = J^ 

4. p = 


j K 

ds 


dcp 2 7^9 


Euler Intrinsic Equation 


5. 

6 . 


— = cos cp = cos(AT-s 2 ) and — = sincp = sin(AT -s 2 ) 

ds ds 

dx coscp j dy . sincp 

— = p coscp = — , — j— and — = p ■ sincp = — . — 

d cp 2V^v9 d s 2v^v9 


Here is a challenge for the DUG-Members: 



In Levies paper ([14] on page 22) I found this pretty figure presenting osculating cubic, 
clothoid and elastica. Can you reproduce this figure? The parameter form of the elastica is 
given below, Josef. 


s s 

J COSCa-SINCb-s)) ds, J SIN(a- SIN(b ■ s)) ds 
_ 0 0 


See: David Halprin , River Meander and Elastica 
in DNL#39. 


Christine Kova’s Students 
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Christine Kova’s Students: 

This is one of the examples we produced in 
a workshop (DERIVE & GeoGebra) several 
years ago. 

One can underlie background pictures of 
ready made stitching (scans from books, 
pictures from the web, or own designs). 

Students learn to apply parameter form and 
families of segments and arcs. 



Student’s works: 
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3D-patterns fixed in acryl cubes: 



This is a more complex pattern. I took it as background picture in DERIVE and you can see 
left below the first step modelling the figure. Josef 







